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XMM- Newton observations of two X-ray bright galaxy 
groups — pushing out to rsoo 
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ABSTRACT 

We report on the results of XMM-Newton observations of the z — 0.18 galaxy group 
WARP J0943. 7+1644 and the nearby poor cluster WARP J0943. 5+1640 at z = 0.256. 
Tracing the X-ray gas out to ~ 425/i _1 and ~ 500/i _1 kpc in the two systems, 
corresponding to roughly 70 per cent of the estimated virial radii, we find the sur- 
face brightness profile of both groups to be well described by standard /3-models 
across the entire radial range, but with a significantly lower value of (3 in the cooler 
WARP J0943. 7+1644. A Navarro-Frenk- White (NFW) model for the gas density can 
also provide a good fit to the surface brightness data, but a large scale radius is 
required, and gas-traces-mass is strongly ruled out. Both systems fall close to the ob- 
served Lx — T relation for clusters. Gas mass fractions increase with radius out to the 
radii probed, reaching extrapolated values at the virial radius of 0.08 and 0.09/i~ 3 / 2 , 
respectively. The latter is concordant with results obtained for more massive systems. 
Gas entropy profiles show evidence for excess entropy out to at least rsoo, and are 
inconsistent with predictions of simple preheating models. This study emphasizes the 
need for observing galaxy groups out to large radii, to securely establish their global 
properties. 
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1 INTRODUCTION 

Groups of galaxies are of fundamental importance to our 
understanding of the mass distribution in the Universe. The 
majority of galaxies are found in groups (Tully 1987), and 
the hot intergalactic gas in groups appears likely to represent 
the dominant component of the Universal baryon content 
(Fukugita, Hogan & Peebles 1998). However, observation of 
this gas presents a problem. The small size of groups, com- 
pared to rich clusters, results in lower X-ray surface bright- 
ness, with the result that whilst cluster gas in some rich 
systems has been mapped out to approximately the virial 
radius r2oo, defined here as the radius containing a mean 
overdensity of 200 with respect to the critical density, the 
emission from galaxy groups is typically lost in the X-ray 
background at r<r2oo/3. 

The seriousness of the problem is aggravated by the 
fact that surface brightness profiles in poor systems are usu- 
ally much flatter than those in clusters (Ponman, Cannon & 
Navarro 1999; Sanderson et al. 2003), probably due to the 
action of non-gravitational processes such as galaxy winds. 
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Adopting the popular /3-model parametrization for the X- 
ray surface brightness S, 

S(r)=So[l + (r/r c ) 2 }- 30+1/2 , (1) 

galaxy groups are generally found to have /3<0.5 (Helsdon 
& Ponman 2000), in comparison to f3 ~ 0.67 in clusters (e.g. 
Jones & Forman 1984; Mohr et al. 1995). Hence both the 
luminosity and gas mass formally diverge at large radius in 
groups. 

The low X-ray surface brightness of galaxy groups has 
so far precluded the detection of the hot gas at large radii 
which should dominate their baryon content. It is important 
to derive observational constraints on the properties of this 
gas to larger radii, for a number of reasons. One can find 
many examples in the literature of authors concluding that 
galaxy groups have low gas mass fractions (defined as the ra- 
tio of gas mass to total mass), on the basis of an X-ray anal- 
ysis extending to only r ~ 200 kpc. While it is known that 
galaxy groups contain large amounts of dark matter, it is not 
yet clear whether their mass-to-light ratios and total baryon 
fractions are indistinguishable from those of clusters, since 
they have not been measured to large radii. Given the domi- 
nant contribution of groups to the mass function of virialised 
systems, this has a substantial impact on attempts to eval- 
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uate the observed mass and baryon content of the Universe. 
It is not clear how properties measured within ~ r2oo/3 
should be extrapolated to 7-200 to obtain integrated values 
for whole systems and provide 'templates' for the interpreta- 
tion of other poor systems. For example, it is not obviously 
appropriate to fit a /3-model and use this to extrapolate the 
X-ray luminosity or gas density distribution. The unphysical 
divergence at large radius discussed above, implies that such 
a profile must turn over at some radius, and one needs to 
know where. In addition, the /3-profile has little theoretical 
justification and better descriptions might be at hand. For 
example, Navarro, Frenk & White (1995; NFW hereafter) 
find from cosmological simulations that out to r2oo, profiles 
of both dark matter and gas density in clusters are better 
fitted by a model of the form 

P{r) K (r/r.)(l+r/r.)» (2) 

(see also Moore et al. 1998). Here r 3 is a characteristic scale 
radius related to the mass of the system. However, these 
simulations do not allow for extra physics, such as galaxy 
winds and radiative cooling, which may radically alter the 
gas profiles, especially in low-mass systems. There is no sub- 
stitute for actually observing how the gas behaves at large 
radii. 

With the arrival of XMM-Newton, we have the oppor- 
tunity for the first time to measure the X-ray properties of 
galaxy groups to large radius. Here we report on the re- 
sults of XMM-Newton observations of the galaxy groups 
WARPJ0943. 7+1644 (hereafter referred to in the text as 
WJ943.7) and WARPJ0943. 5+1640 (hereafter WJ943.5 in 
the text), aimed at mapping and characterizing the X-ray 
gas out to a considerable fraction of 7-200 • These X-ray 
sources are close on the sky and were targeted in a single 
XMM pointing; they were both detected serendipitously in 
a 9 ks ROSAT PSPC pointing as part of the WARPS X-ray 
cluster survey (Scharf et al. 1997; Jones et al. 1998). The 
extended X-ray emission of the systems was also confirmed 
by their subsequent detection in the CfA 160 deg 2 survey 
(Vikhlinin et al. 1998; Mullis et al. 2003) and in the RIXOS 
survey (Mason et al. 2000). In the optical follow-up of the 
former, spectroscopic redshifts of 0.180 and 0.256 were de- 
rived for WJ943.7 and WJ943.5, respectively. 

In iJ5]we describe the data preparation, ^ deals with 
the methods and results of the spatial and spectral analysis 
of the groups, while [|H describes results and implications 
of the mass analysis. A summary is presented in 

We adopt a cosmology with Q, m — 0.3 and J1a = 0.7, 
and write the Hubble parameter as Ho = 100ft km s _1 
Mpc" 1 . 



2 OBSERVATIONS AND DATA REDUCTION 

WARPJ0943. 7+1644 and WARP J0943. 5+1640 were ob- 
served by XMM-Newton in two separate observations for 
a total exposure time of 22 ks. All detectors were in Full 
Frame mode, and a medium optical blocking filter was used 
due to the presence of nearby bright stars. 

Data preparation was carried out using the XMMSAS 
v5.3 software. Screening for periods of high background was 
performed by inspection of lightcurves extracted in the 10- 



Figure 1. Combined (pn+2MOS) 0.4-2.5 keV photon count im- 
age after initial data filtering. The central 15 X 15 arcmin 2 is 
shown, smoothed with a Gaussian of 10 arcsec (FWHM). North is 
up and East is to the left. WJ943.7 is the central extended source, 
while WJ943.5 is seen roughly 5 arcmin to the South-west. 

12 keV (for MOS) and 12-15 keV (pn) bands. Periods with 
count rates in excess of 0.15 cts s _1 for MOS1/2 and 0.22 
cts s _1 for pn were removed. The second of the two ob- 
servations was heavily affected by periods of high particle 
background. Screening for this left 7.7, 14.5, and 14.4 ks of 
useful exposure time for the pn, MOS1, and MOS2 cameras, 
respectively, in the first observation, while the correspond- 
ing values for the second observation were 1.0, 4.7, and 4.6 
ks. 

2.1 Imaging 

Fig. Q shows an image in the 0.4-2.5 keV band after the 
initial data filtering. Bright point sources were subsequently 
identified with XMMSAS and excised using an exclusion ra- 
dius of 25 arcsec. Images were background subtracted prior 
to analysis as follows. For each instrument and observation, 
a blank sky background event file (Lumb 2002) was gener- 
ated which had undergone identical lightcurve filtering to 
the data, and had been coordinate transformed to match 
the aspect solution of the latter. This ensured that the same 
CCD regions were considered when comparing background 
levels in our data and the blank sky files. The blank sky 
background was then determined in the 9-12 keV band 
within an annulus of inner radius 9 arcmin and outer ra- 
dius 11 arcmin, and an analogous estimate was made for 
our data. Comparing the two background levels yielded a 
scaling factor / between the high-energy background in the 
current observations and in the blank sky data. The total 
background of the latter was then scaled accordingly and 
subtracted from the data. 

At energies E< 1.5 keV, the resulting images may still 
contain a contribution from the soft X-ray background, 
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which is known to vary considerably with sky position and 
is not necessarily well represented by the blank-sky data for 
a given observation. The soft (E < 2 keV) photon excess of 
the resulting images relative to the blank-sky background 
was therefore determined in a large-radius (9-11 arcmin) 
'source-free' annulus and subtracted from the images after 
correction for vignetting. This excess was found to be neg- 
ative, i.e. a deficit of soft emission relative to the blank sky 
files is present in the data, constituting on average ~11 per 
cent of the total background over the entire field of view (a 
deficit might be expected, as the value of absorbing hydro- 
gen column density Nh is on average a factor of ~2 lower 
in the Lumb background data than at the position of our 
pointing) . 

Note that all subsequent analysis, including investiga- 
tion of source spectra, was restricted to energies E > 0.4 
keV, for which the effective-area difference is negligible be- 
tween the thin optical filter used for the Lumb (2002) data 
and the medium filter applied in our observations. 



2.2 Spectra 

A similar approach to the above was used for the removal of 
background in extracted spectra. Using the scaling factors 
/ derived above, source spectra were created, and scaled 
blank-sky spectra extracted from the same (detector) re- 
gion were subtracted for each instrument and observation. 
Large-radius (9-11 arcmin) source-free spectra were also 
created from the data and associated scaled blank-sky spec- 
tra were subtracted, the result reflecting the spectrum of any 
remaining soft background component due to the particu- 
lar position on the sky. These residual spectra were finally 
also subtracted from the source spectra. All results from 
spectral fitting presented in this paper are based on this 
method for background subtraction, using the Lumb (2002) 
data as background templates. However, as a check on the 
robustness of our results, comparisons were made with re- 
sults from two other methods for background subtraction. 
One is identical to the above, but uses another set of blank- 
sky background data also assembled from pointed observa- 
tions (Read & Ponman 2003). These background data have 
the advantage of having been obtained with the medium 
blocking filter, allowing a test for any effect which might 
arise from using the thin-filter Lumb (2002) data for back- 
ground subtraction in our medium-filter observations. The 
other method involves a simple local background subtrac- 
tion, extracting background spectra in a 9-11 arcmin an- 
nulus in the source data. This ensures that the effect of 
any discrepancy between absorbing hydrogen column den- 
sity and optical filter in source and background data can be 
tested, but it does not account for chip-to-chip variations in 
response or background level. We find that these two alter- 
native methods produce results for the derived parameters 
which are consistent within la with those obtained from the 
adopted approach. We therefore take our spectral results to 
be fairly robust. 



3 DATA ANALYSIS 

3.1 X-ray/optical identification and morphology 

As is readily apparent even in the raw imaging data (Fig.0, 
there is evidence for a third extended extended source in the 
field in addition to the two WARPS groups: To the South- 
east of the optical axis, faint extended emission is seen which 
was apparently not detected in R OS AT pointings and is con- 
sequently not included in any ROSAT-based source cata- 
logues. We shall refer to this source as XMMJ0943. 9+1641. 

To determine the detailed X-ray morphology of these 
three systems as seen by XMM, images were produced for 
each instrument and observation in the 0.4-2.5 keV band 
(chosen to optimize S/N on the basis of extracted source 
spectra, see § 13.31 . from which a combined, background 
subtracted, and exposure corrected image was made. Fig. 
shows resulting X-ray surface brightness contours of the gen- 
eral region surrounding WJ943.7 and of the three extended 
sources in the region. Individual source contours are overlaid 
on optical images; where possible, we have used the follow- 
up imaging data from the WARPS survey itself. 

To aid in investigating the nature of optical sources in 
the fields, and since spectroscopic results for the two groups 
are not available in the literature, the online near-infrared 
(JHK bands) catalogue of the 2 Micron All Sky Survey 
(2MASS) and its listing of associated optical sources in the 
Sloan Digital Sky Survey was employed. ^From the aper- 
ture magnitudes of the sources, colour-magnitude diagrams 
for all sources within 0.7ft -1 Mpc of the X-ray centroid were 
produced (approximately the virial radii of the systems, see 
corresponding to angular scales of 5.5 arcmin at z = 0.18 
and 4.1 arcmin at z = 0.256. For XMMJ0943.9+1641 a ra- 
dius of 4.0 arcmin was conservatively used). Colours were 
compared to those expected for a population of passively 
evolving early-type cluster galaxies of various metallicities 
at the appropriate redshift (Kodama & Arimoto 1997). A 
formation redshift of Zf = 3 was assumed, but results are 
essentially insensitive to realistic values of this parameter. It 
is worth keeping in mind when distinguishing between point- 
like and extended NIR sources in the 2MASS catalogue that 
the point spread function (PSF) of the 2MASS telescopes is 
2.5 arcsec at FWHM, corresponding to ~5/i _1 and ~7/i _1 
kpc at the redshifts of WJ943.7 and WJ943.5, respectively. 
At faint NIR fluxes, a source classified as pointlike could 
thus potentially be a galaxy, if only the central bulge is de- 
tectable. 



3.1.1 WJ943.7 

WJ943.7 was selected for this study in part on the basis of 
its seemingly relaxed morphology in ROSAT PSPC observa- 
tions. The overall structure of the X-ray emission is indeed 
seen from Fig. [5] to be fairly symmetric, supporting the as- 
sumption that the system is reasonably relaxed. 

The position of the X-ray peak and centroid at (a — 
09 h 43 m 44 s . 6, S = +16°44'25") coincides with that of an 
optical source, clearly identified as a galaxy upon inspec- 
tion of the image. No optical identification of this galaxy 
could be secured from the LEDA or NED online databases. 
Its ii-band apparent magnitude in the 2MASS catalogue is 
rriR — 16.0, making it the brightest galaxy in the field shown 
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Figure 2. Top left: 0.4—2.5 keV surface brightness contours of the combined (pn+2MOS) background-subtracted and exposure-corrected 
image, binned to a spatial scale of 4 arcsec pixel -1 and smoothed with a 20 arcsec Gaussian (FWHM). Dashed circles mark the exclusion 
regions invoked for the surface brightness and spectral analysis of WJ943.7 (see text). Top right: Smoothed and linearly spaced contours 
(0.6-3 Xl0~ 5 cts s" 1 pixel" 1 ) of XMMJ0943.9+1641 overlayed on a 5 X 5 arcmin 2 optical image from the Digitized Sky Survey. Bottom 
panel: Contours for WJ943.7 (0.6—6 in the above units) and WJ943.5 (0.8—8), overlayed on WARPS survey /-band and H-band image, 
respectively, taken with the 1.3-m McGraw-Hill telescope at MDM Observatory (images courtesy of G. Wegner and H. Ebeling). 



in Fig. |21 Assuming an early-type galaxy spectrum and ac- 
counting for k— correction (Poggianti 1997) and Galactic ex- 
tinction (Schlegel, Finkbeiner & Davis 1998), this translates 
into My ~ —22.9 (h = 0.7), which is indeed consistent with 
expectations for a bright elliptical. The source displays ex- 
tended (7 arcsec) NIR emission in the 2MASS data and we 
find that its B-K, J-K, and H-K colours are perfectly consis- 
tent with an elliptical galaxy at the spectroscopic redshift. 
A total of four additional extended NIR sources are found 
within 0.7/i _1 Mpc, of which two have colours very similar 
to those of the central source. These four sources span an 
i?-band apparent magnitude range of ma = 16.8 — 18.1. The 
bright source to the North-west of the X-ray centroid is a 
star. 



3.1.2 WJ943.5 

This source was similarly detected in both the CfA 160 deg 2 
(Vikhlinin et al. 1998) and RIXOS (Mason et al. 2000) sur- 
veys. In contrast to WJ943.7, the X-ray emission of WJ943.5 
clearly seems elongated in the NE-SW direction. The ques- 
tion is whether this apparent elongation is partly or entirely 
due to the presence of two nearby point sources (most clearly 
seen in the top left of Fig. |5J , in particular the one East of 
the emission peak. We will return to this issue in H3.2l and 
H3.3I Nevertheless, the overall structure of emission appears 
slightly less relaxed than for WJ943.7. 

Again, the position of the X-ray peak at (a = 
Q9 h 43 m 31".0,S = +16°39'57") (displaced ~ 15 arcsec to the 
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W with respect to the X-ray centroid) coincides with an 
optical source. The deep optical exposure shown in Fig. [2] 
allows us to readily identify it as an early-type galaxy. It is 
not listed in LEDA or NED but features an extended coun- 
terpart in the 2MASS data, this being however the only 
extended NIR source within 0.7/i _1 Mpc. Its colours are in 
good agreement with those expected for a cluster elliptical 
at z = 0.2 ~ 0.3. Investigation of the image picks up at least 
a handful of more galaxies within the region marked by the 
X-ray contours, none of which have a NIR counterpart in 
the 2MASS data. This is not surprising, however, as the 
2MASS nominal survey completeness limits are mj = 15.8, 
ran = 15.1, and ulk = 14.3 mag (for pointlike sources), 
thus rendering the detection of L < L* galaxies at z = 0.256 
rather unlikely. The two optically bright sources within ~ 1 
arcmin East of the centre are stars. 

Regarding the two X-ray point sources, faint optical 
sources are present at the position of both, but neither dis- 
play detectable NIR emission in the 2MASS data. As will be 
discussed, X-ray spatial and spectral analyses suggest that 
the Western source is truly pointlike, whereas the Eastern 
source might contain an extended component. In all subse- 
quent analysis both sources were masked out, adopting ex- 
clusion radii of 25 arcsec and 40 arcsec for the Western and 
Eastern source, respectively. These choices will be justified 
in fl3~2landl3~31 



3.1.3 XMMJ0943. 9+1641 

As already mentioned, this source was not identified in 
ROSAT pointings. We regard the identification of the source 
in these data as secure, with the innermost contour in Fig. [5] 
being 8a above the residual background. While no clear 
emission peak is identifiable, the X-ray centroid is found to 
be at (a = 09 h 43 m 58 ! :5, 5 = +16°41'17")- The structure of 
emission seems quite disturbed, but deeper observations are 
required to firmly establish its X-ray morphological details. 

Inspection of the 2MASS catalogue and Fig. [2] revealed 
evidence for four NIR/optical sources within the region cov- 
ered by the displayed X-ray contours, of which only one 
is classified as extended in the infrared (roughly 0.7 arcmin 
from the X-ray centroid). None of these sources stand out as 
a bright central galaxy. They all have very similar colours, 
consistent with those of a cluster population at redshifts 
between 0.1-0.4 (with perhaps the upper value being more 
likely, given the optical and X-ray faintness of the region rel- 
ative to the two WARPS groups). A second extended NIR 
source is present, but at a centroid distance of ~ 3 arcmin 
to the NW, it might rather be associated with WJ943.7. 



3.2 Surface brightness distribution 

3.2.1 One- dimensional profiles 

The radial surface brightness profiles of the two WARPS 
groups were initially parametrized by fitting standard /3- 
models, eq. Q, to the azimuthally averaged exposure- 
corrected emission, including a spatially uniform residual 
background. To minimize contamination from the other two 
extended sources in the field, exclusion regions for each of 
these were defined at their respective emission centroids (see 



Fig.|5|for an example). A residual background level was de- 
termined, and fitting was done through adaptively binning 
the brightness profile such that each radial bin contained 
emission at >4.5cr above this background; this allowed an 
estimate of the maximal extent of emission while at the same 
time ensuring the presence of at least 20 source counts per 
bin (justifying the use of \ 2 statistics for fitting). 

The resulting radial profiles and best-fitting /3-models 
are shown in Fig. [3] and in Table^below we list the derived 
best-fitting parameters for both groups. Emission at >4.5<r 
is detected out to at least 3 arcmin in both cases, and the fits 
are acceptable, with reduced x 2 = 0.94 and 0.86 for 8 and 
10 degrees of freedom (d.o.f.), respectively. Neither system 
shows evidence for any significant large-radius departure 
from the /3-model or excess emission in the centre indicating 
a cooling flow. Note that WJ943.5 is detected to a smaller 
radius than WJ943.7, although it displays a higher central 
surface brightness. This reflects its larger /3 and hence fainter 
fluxes at sufficiently large radii. For XMMJ0943.9, emission 
at > 4.5(7 was detected out to ~ 75 arcsec but statistics were 
too poor to obtain useful constraints from a /3-model fit. 

The gas density profiles in groups must break from a 
power law decline at sufficiently large radii. Cosmological 
simulations which include only gravity and shock heating, 
produce clusters in which the gas density distribution follows 
that of the dark matter, and steepens progressively towards 
an r~ 3 dependence at large radii (NFW; Eke, Navarro & 
Frenk 1998; Frenk et al. 1999). Motivated by this, we also 
fitted a gas density profile of the NFW form, eq. 0, to the 
surface brightness profiles. Since there is little evidence for 
non-isothermality in the two groups ( H3.3L we assumed X- 
ray emissivity tx n\ , where n e is the electron number density. 
For a profile of the form of eq. (2) , the resulting line of sight 
projection Sf^p^Y is not analytically tractable. Hence the 
fit was done by performing x 2_mm i m i za tion on a grid of 
numerical profiles. 

At large radii (r S> r s ), Snfw oc r -4 , which is clearly 
much steeper than the dependence revealed by Fig. El Large 
scale radii r s are therefore anticipated in the fits, as are 
large discrepancies in the inner regions where the observed 
flat density cores are incompatible with the steep NFW be- 
haviour. Consequently, the data points inside the best-fitting 
core radii of the /3-models were excluded in the fits. For 
both groups, the isothermal NFW fits are good, yielding 
fit qualities comparable to those obtained using /3-models, 
the fit being slightly better for WJ943.5 and slightly worse 
for WJ943.7 than a /3-model fit. Thus, assuming isothermal 
gas, the data are in accordance with models where the gas 
density shows a progressively steepening power-law index, 
as it does for an NFW model. As will be shown in i j4.ll 
the X-ray emissivity of the groups is independent of tem- 
perature to within ~ 5 per cent for the relevant range of 
parameters involved, so the above result applies even if the 
gas is moderately non-isothermal. However, the scale radii 
obtained are in all cases large: for WJ943.5 r s ~ 630h~ kpc 
and for WJ943.7 r s ~ 850/i _1 kpc (even larger than the de- 
rived r2oo, cf. §EJ. This corresponds to small concentration 
parameters c = r2oo/r s < 1, and is in apparent conflict with 
expectations from numerical simulations (NFW; Eke et al. 
1998; Moore et al. 1998; Lewis et al. 2000) and other ob- 
servations (e.g. Pratt & Arnaud 2003) which typically find 
r s ~ 0.1 — 0.25r2oo (i.e. c w 4—10) for the dark matter distri- 
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Figure 3. Top: Radial surface brightness profiles of WJ943.7 (left) and WJ943.5 (right), binned such that each bin contains emission at 
greater than 4.5cr significance. Also shown is the best-fitting /3-model in each case. Error bars are la. Bottom: Corresponding confidence 
contours (la, 2a, 3a) in the (f), r c ) plane. Solid (dark) contours describe those pertaining to the group itself, whereas dashed (light) 
contours describe those of the other group, overlaid for comparison. The same scales have been used in each panel for ease of comparison. 



bution. Given that we are dealing here with the gas density, 
this most likely simply reflects the fact that the gas distri- 
bution is much less centrally concentrated than that of the 
dark matter, as will be substantiated in Models in which 
gas traces total mass are hence clearly ruled out for both 
our groups. This is also supported by the results of Ettori 
& Fabian (1999), who use a similar approach for a sample 
of massive clusters, and find a mean ratio between NFW 
gas scale radius r s and /3-model core radius r c of r s /r c ~ 3, 
whereas we find r a /r c >15 for WJ943.7 and r s /r c >5 for 
WJ943.5. Importantly, the large scale radii resulting from 
the fits also show that the derived surface brightness pro- 
files are inconsistent with a steepening in the logarithmic 
density slope r\ = — (dlnn e /dlnr) towards r\ ~ 3 at large 
radii. In fact, the best-fitting NFW models produce slopes 
of n = 1.66 and 1.80 for WJ943.7 and n = 1.87 and 2.08 for 
WJ943.5 at the derived extent of emission r ex t (see Tabled 
below) and at r2oo (3D, respectively. For both groups, the 
concluding point is therefore that isothermal NFW models 



only provide a good description of the data if the associ- 
ated scale radii are so large as to make the resulting surface 
brightness profiles rather similar to /3-models. 

3.2.2 Two-dimensional profiles 

Although the X-ray emission of both WARPS groups is rea- 
sonably symmetric and is well fitted by a one-dimensional 
/3-profile, two-dimensional /3-models were also fitted to the 
emission, using SHERPA, in order to further investigate the 
surface brightness structure. This approach has the advan- 
tage that the effects of the asymmetric EPIC PSF can be ac- 
curately taken into account in the model fit. However, whilst 
2-D fitting takes better advantage of the available spatial in- 
formation, the image counts do not satisfy the requirement 
of Gaussian errors which underlies the \ 2 statistic (the num- 
ber of counts in some data bins may be very low). Following 
Helsdon & Ponman (2000), the maximum-likelihood based 
Cash statistic (Cash 1979) was used instead. Vignetting was 
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Table 1. Best-fitting values from 1-D and 2-D surface brightness fitting (0.4-2.5 keV). All uncertainties are la for one interesting 
parameter. Physical scales are computed for distances corresponding to the spectroscopic redshifts of the groups. Radial extents r ex t 
specify the outer radius of the outermost radial bin containing a 4.5<r detection. Position angles are given counter-clockwise from the 
North. The goodness of fit for 2-D models is as described in the text. 
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taken directly into account in the fit through including the 
exposure map in the modelling of the source. 

The 2-D fit results are also presented in Table For 
each source, the PSF model employed for convolution with 
the source model was produced at the relevant detector po- 
sition at a photon energy of 1 keV using XMMSAS. Additional 
free parameters compared to the 1-D fit were the eccentric- 
ity e = (1 — & 2 /a 2 ) 1//2 , position angle 8, and x- and j/-centre 
of the X-ray gas distribution. 

No estimator of the goodness of fit is immediately avail- 
able when using the Cash statistic. Hence to address the is- 
sue of fit quality, we adopted a Monte Carlo approach, gen- 
erating for each source 1000 artificial images based on the 
best-fitting model. In each case, Poisson noise was added to 
the input model image, which was subsequently convolved 
with the PSF model and vignetted. Model fits were then per- 
formed to these images in the same manner as for the real 
data, recording for each fit the value of the Cash statistic, 
and finally fitting a Gaussian to the resulting distribution of 
these values. Comparing the Cash statistic, C, derived for 
the real data to this distribution, a measure of the goodness 
of fit was obtained as the number of standard deviations a 
separating C from the centre of the distribution. This value 
is also presented in Table Q and it is seen that both fits are 
reasonably good, being less than la from the distribution 
centre. Note also that confidence intervals can be assessed 
as for x 2 statistics, since differences in the Cash statistic are 
X 2 -distributed. 

For both groups, the 1-D and 2-D results agree well 
within the errors. For WJ943.7, the derived /3-value is typ- 
ical for groups and poor clusters with X-ray temperatures 
below 1-2 keV. The X-ray morphology is consistent with 
zero eccentricity, and 9 is therefore unconstrained for this 
source. WJ943.5, in contrast, appears to display a larger 
X-ray core and a /3-value more typical of richer clusters 
(T>3 keV). The derived eccentricity of its X-ray emission, 
e ~ 0.2, is consistent with typical values found for other 
groups and clusters (Mohr et al. 1995; Helsdon & Ponman 
2000). R OS A T results for other groups (Mulchaey & Zablud- 
off 1998) have shown that the optical light of the brightest 



group/cluster galaxy is often well aligned with the position 
angle of the X-ray emission. In the case of WJ943.5, how- 
ever, the X-ray position angle 9x ~ 80° seems not particu- 
larly well aligned with the .R-band optical light of the cen- 
tral galaxy, which shows 8 op t ~ 50°. In the central regions, 
within 1 arcmin of the emission peak, the X-ray contours are 
nearly circular (e ~ 0.02), but the best-fitting position angle 
of 9x ~ 20° is effectively unconstrained and so is consistent 
with both of the above values. 

The difference in r c , and particularly /3, seen for the 
two groups is interesting, inasmuch as the systems have rea- 
sonably similar temperatures ( H3.31 and hence masses. To 
investigate the significance of the difference, we also plot in 
Fig- 01 the derived confidence contours in the {/3,r c ) plane 
for both groups. Although the la contours marginally over- 
lap, it is clear that the two profiles are not consistent at la. 
For WJ943.5, any combination of r c and (3 leading to a typ- 
ical group value of (5 < 0.5 is only marginally allowed within 
3cr confidence. It could of course be conjectured that some 
process has lead to a significantly larger X-ray core for this 
group, which, given the quality of the data, naturally forces 
/3 towards high values in the fits. Another possibility is re- 
lated to the fact that when masking out the Eastern point 
source seen in Fig. to the standard radius of 25 arcsec, 
resulting residuals from 2-D surface brightness fitting sug- 
gest the source is not quite pointlike (whereas the W source 
is well accounted for by point source emission). This was 
confirmed by explicitly including, in the overall fit to the 
group emission, a point source model at the position of the 
E source, the result still showing systematic residuals around 
the source. The results listed in Tableware therefore based 
on an exclusion radius of 40 arcsec for this source, estimated 
from the 2-D fitting residuals to mask out the excess emis- 
sion to sufficient accuracy. Furthermore, as will be described 
below, spectral results also suggest that the Eastern source 
is not pointlike. If indeed a 'separate' extended component 
is present at ~ 1 arcmin East of the group emission centroid, 
this might also drive r c , and hence (3, towards high values 
in the fits. 

We note that an eccentricity of ~ 0.2 persists when 
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explicitly including the Eastern source in the 2-D modelling, 
either as pointlike or extended. This strongly suggests that 
the derived elongation and position angle of the group X- 
ray contours are intrinsic and not caused by excess emission 
from an extended source embedded in the group emission. 

3.3 X-ray temperatures 

For the spectral analysis, integrated spectra of the WARPS 
groups out to r ext were extracted for each instrument and 
observation, accumulating the spectra in bins containing 
at least 20 counts. For each source, all six background- 
subtracted spectra were then fitted simultaneously in the 
0.4-4 keV band, using an absorbed mekal model in xspec 
vll.l, with the redshift fixed at the spectroscopic value and 
Nh fixed at the Galactic value of 3.2 x 10 20 cm -2 (Stark et al. 
1992). Free parameters were thus temperature, metal abun- 
dance, and spectral normalization. Results are presented in 
Tabled with metallicities derived using the solar abundance 
table of Anders & Grevesse (1989). Acceptable fits were ob- 
tained, with best-fitting temperatures and abundances char- 
acteristic of X-ray bright groups or poor clusters. 

For WJ943.5 it was found that the spectrum could only 
be adequately fitted by one thermal component, provided 
the Eastern point source was excised out to at least 40 arc- 
sec, once again suggesting that this source is extended rather 
than pointlike. The combined spectrum of the two 'point' 
sources seen in the group is nevertheless well modelled by 
a power- law of slope 1.8 ± 0.2 added on top of the ther- 
mal plasma spectrum of the intragroup medium. However, 
masking out the Western source and modelling the Eastern 
alone along with a thermal plasma representing the group 
emission, we find that in all cases marginally better fits are 
obtained when replacing the power-law model with another 
mekal thermal plasma at the group redshift, fixed at a metal- 
licity of O.3Z0 . We find no evidence that the source is rather 
described by multiple power-laws or a single power-law with 
internal absorption. Allowing the redshift of the mekal model 
to vary yields z in the interval 0.2-0.3, i.e. consistent with 
the group redshift. These conclusions apply independently 
of the exact choice of fitting region, and seem to suggest 
the presence of a separate blob of plasma at a temperature 
roughly 0.5 keV hotter than the surrounding medium. Due 
to the limits on the data, a hardness map does not provide 
any additional information in this context. 

For XMMJ0943.9, we also attempted to fit a thermal 
plasma model to the 0-75 arcsec spectrum, keeping temper- 
ature, redshift, and normalization as free parameters, while 
fixing the abundance at O.3Z0. A mekal model provides a 
good fit (x 2 = 20.9 for 23 d.o.f.), yielding T = 1.40+°; 2 g (lcr 
errors) and a redshift of z = 0.33^0 07 (90 P er cent confi- 
dence). With both parameters reasonably well determined, 
this suggests that XMMJ0943.9 is a yet more distant galaxy 
group. The derived redshift is in good agreement with that 
expected from the optical-NIR colours of sources in this re- 
gion, cf. H3.ll We derive an unabsorbed 0.5-2 keV flux of 
~ 3 x 10~ 14 ergs cm -2 s _1 , corresponding to a luminosity of 
1.0 x 10 43 h~ 2 ergs s _1 at the nominal best-fitting redshift. 

Due to the limited statistics available (we record 
~1500 and ~1100 background-subtracted source counts 
for WJ943.7 and WJ943.5, respectively, within the non- 
excluded regions), a detailed study of temperature and 



Table 2. Results of mekal fits to the integrated 0.4—4 keV spec- 
tra, with Nh fixed at the Galactic value of 3.2 X 10 20 cm -2 . 
Luminosities are given in the 0.5-2 keV band in units of 10 43 h~ 2 
ergs s _1 (see text for details). Errors are 1<t for one interesting 
parameter. 
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abundance gradients in the intra-cluster medium (ICM) is 
not feasible. Although a hardness map does not reveal ev- 
idence for any significant temperature variations within ei- 
ther group, an examination of the gross structure of the ICM 
temperature distribution T(r) in the two groups was never- 
theless attempted by extracting spectra in two or three con- 
centric annuli, containing roughly equal numbers of source 
counts. The spectra were then fitted as above, with abun- 
dances fixed at the value derived from the fits to the inte- 
grated spectra inside r ext (consistent results for both T and 
Z were obtained when the latter was allowed to vary). In 
all bins good fits are obtained, and results for two bins are 
shown in Fig. [I] Although both groups show tentative evi- 
dence for a mild (~0.5 — 1 keV) temperature decline with 
radius, the ICM is consistent with being isothermal within 
the 90 per cent confidence (1.65<r) errors in each bin. Specif- 
ically, there is no indication of the presence of cool cores in 
the central regions. This is consistent with the lack of ex- 
cess emission in the group cores relative to the best-fitting 
isothermal /3-model. When spectral properties are derived 
using three rather than two annuli, these are, of course, sub- 
ject to larger statistical uncertainties, but are again consis- 
tent with isothermality, and show no central cooling. 

3.4 X-ray fluxes and luminosities 

From the best-fitting mekal models for the integrated spec- 
tra, resulting 0.5-2 keV unabsorbed fluxes for WJ943.7 
(WJ943.5) are 0.96l°;^(1.12l^ 4 ) x 10" 13 ergs cm" 2 s" 1 
within the regions adopted for spectral fitting. Errors were 
determined from the fractional errors on the spectral nor- 
malisation. Correcting the WJ943.5 flux inside the spectral 
region for removed area, assuming the nominal best-fitting 
2-D /3-model in Table[T| the resulting flux is 1.43±o;i| x 10~ 13 
ergs cm -2 s , where the quoted la errors are from the spec- 
tral flux determination alone. 

At the relevant redshifts, the derived fluxes within r ex t 
correspond to unabsorbed 0.5-2 keV restframe luminosi- 
ties of L x = 0.57l°;° 7 x 10 43 /i" 2 ergs s" 1 (WJ943.7) and 
1.97^ 25 x 1q43 h ^ er S s s_1 (WJ943.5, corrected for re- 
moved area as above) , as listed in Table |21 In Fig. |5] the 
best-fitting Tx and the bolometric (0.01-15 keV rest-frame) 
luminosities within the radii adopted for spectral fitting are 
compared to those of the low-redshift (z < 0.04) T > 1 keV 
groups in the samples of Helsdon & Ponman (2000) and 
Xue & Wu (2000). Note that both depicted L x - T rela- 
tions for groups {Lx oc r 4 9±0 8 , Helsdon & Ponman 2000; 
L x oc T 5 57±179 , Xue & Wu 2000) have been derived using 
groups at 'all' T, whereas only groups with T > 1 keV are 
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Figure 4. Projected radial variation of the X-ray gas temperature in WJ943.7 (left) and WJ943.5 (right). Crosses mark the best-fitting 
values in each bin, with surrounding solid (dotted) lines showing the 90 per cent (ltr) confidence intervals. Dashed horizontal lines 
show the best-fitting value from the integrated spectrum, with the associated la errors marked by dotted lines. The dot-dashed curves 
represent the polytropic models discussed in the text. 



included in Fig. |S] Note also that while the groups in the 
Helsdon & Ponman sample were analysed in a coherent man- 
ner, data for the Xue and Wu groups were compiled from a 
number of different sources. In neither case was the X-ray 
emission of the groups derived within a consistent radius, 
be it physical or as a fraction of r2oo, and thus a rigorous 
comparison to our results cannot be made. Nevertheless, it 
is found that both WARPS groups have luminosities con- 
sistent with those derived for other poor clusters, agreeing 
well with the cluster Lx — T relation of slope 2.79 ± 0.08 
derived by Xue & Wu (2000). This L x -T relation is based 
on a sample of 274 clusters drawn from the literature, the 
large majority of which are at z < 0.5, and agrees well at 
the relevant temperatures (1^T<3 keV), with, e.g., the 
low-to-moderate redshift (z < 0.2) relation derived for the 
T > 1 keV WARPS clusters (Fairley et al. 2000) and with 
estimates for z < 0.1 clusters from BeppoSAX data (Ettori, 
De Grandi & Molendi 2002). 
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4 MASS ANALYSIS 

4.1 Effects of non-isothermality 

It is clear from Fig. Q]that the temperature profiles are not 
resolved with sufficient accuracy to immediately allow de- 
tailed profiles to be made of temperature-based quantities 
such as mass or gas mass fraction. Fig. 2] and the ensu- 
ing discussion shows, however, that there is no evidence for 
significant departures from isothermality in either group at 
any radius probed, and an isothermal distribution is, in- 
deed, the flattest reasonable temperature distribution with 
which the data are consistent. One might argue that T in 
the innermost bin could have been underestimated due to 
the presence of a cooling flow; this would produce a steeper 
temperature profile outside the cooling core and hence a 
larger deviation from isothermality than is immediately ap- 
parent in Fig. [I] While the expected central cooling times 



Figure 5. The derived bolometric luminosities and temperatures 
of the two WARPS groups (diamonds), shown together with the 
values for the T > 1 keV groups of Xue & Wu (2000; triangles) 
and Helsdon & Ponman (2000; squares). Also shown is the nom- 
inal best-fitting Lx — T relation for all groups in the samples 
of Xue & Wu (dashed line), Helsdon & Ponman (dotted line), 
and for the cluster sample of Xue & Wu (solid line). The figure 
assumes an Einstein-de Sitter cosmology with h = 0.5. 

(using the cooling curves of Sutherland & Dopita 1993, for 
a Z = 0.3Z Q plasma) of ~ 7 and ~ 11 Gyr for WJ943.7 
and WJ943.5, respectively, may suggest the presence of 
weak cooling flows, this is not supported by the derived 
surface brightness profiles. Moreover, a number of recent 
studies of low-temperature systems with well-defined tem- 
perature profiles have demonstrated that, outside a possible 
cool core, isothermality may indeed be a good assumption 
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(e.g. Mushotzky et al. 2003; Pratt & Arnaud 2003). In the 
following sections we will therefore assume the ICM to be 
isothermal at the mean temperature derived inside r ex t- It 
is nevertheless instructive to consider the impact of assum- 
ing isothermality, in order to assess the potential error in 
doing so, and in all subsequent analysis we will attempt to 
quantify potential effects of this assumption on our results. 

On the basis of the temperature data of Fig. 2] one 
can assume model profiles for T(r) which (i) are consistent 
with the data, (ii) can be justified with reference to other 
observed and simulated systems, and (iii) can provide an in- 
dication of the systematic errors associated with making a 
single isothermal assumption for T(r). Hence, in the absence 
of detailed temperature data, such model profiles are used in 
the following in conjunction with the isothermal assumption. 
We readily acknowledge the fact that the resulting radial 
variation of various quantities derived from T(r) is, strictly 
speaking, modelled rather than observed, providing an ap- 
proximation to the true profile which is only reliable to the 
extent to which the ICM temperature behaves as assumed. 
As will be seen, the overall conclusions are nevertheless not 
sensitive to specific choices for T(r) allowed by Fig. [I] 

The detailed shape of the temperature profile in our 
groups is largely unconstrained. There is, however, com- 
pelling evidence from both observations and simulations of 
groups and clusters that a polytropic equation of state, 

Tcx^" 1 , (3) 

provides a good description of any radial decline in gas tem- 
perature outside a possible cool core. This is the case both 
for simulated clusters (Lewis et al. 2000; Loken et al. 2002), 
simulated groups/cool clusters (Ascasibar et al. 2003), 
and for observed clusters with well-determined tempera- 
ture profiles (e.g. Markevitch et al. 1998; Ettori & Fabian 
1999; Markevitch et al. 1999; Pratt, Arnaud & Aghanim 
2001; Pratt & Arnaud 2002) or surface brightness profiles 
(Miralda-Escude & Babul 1995). It is clearly also a better 
description than the isothermal assumption in some clus- 
ters with less well-determined temperature profiles (see e.g. 
Sarazin, Wise & Markevitch 1998). Eq. (J^J is consequently 
employed here as a tool to investigate the systematic errors 
associated with assuming isothermality. For the tempera- 
ture data shown in Fig. 2] polytropic models intercepting 
the nominal best-fitting temperature in each bin have in- 
dices of 7 = 1.15 for WJ943.7 and 7 = 1.22 for WJ943.5, 
very similar to values found for observed (Markevitch et al. 
1999; Ettori & Fabian 1999) and simulated (Lewis et al. 
2000; Ascasibar et al. 2003) clusters. Here we shall not at- 
tempt to include errors on these indices; given the available 
photon statistics, they would not be well constrained for 
the present data and their errors would propagate into large 
uncertainties on the resulting total masses. 

If indeed the gas is polytropic, i.e. satisfying eq. @, 
and if p follows a /3-model, the true gas temperature Tt rtle 
will differ from the projected (i.e. measured) one T pro j only 
by a constant factor, 

Ttrue = r[|^(i + 7 )]r(3/3-|) 

T vroj r[f/3(l+7)-|]r(3/3) 

(Markevitch et al. 1999), which amounts to 1.07-1.08 for 
the systems studied here. Strictly speaking, eq. @ is only 
valid assuming emissivity oc n^T a , but in the 0.4-4 keV 



band considered for spectral fits, and over the relevant tem- 
perature ranges (1.3-2.0 keV for WJ943.7 and 1.6-2.9 keV 
for WJ943.5; see Fig. 0J, the emissivity is temperature- 
independent to within 5 per cent and is oc T a to within 
3 per cent for a Z — 0.3ZQ mekal plasma. 



4.2 Total masses and virial radii 

Motivated by the regularity of the X-ray emission in the 
groups, we assume spherical symmetry and hydrostatic equi- 
librium of the X-ray emitting gas. The total gravitating mass 
inside a radius r is then given by 

M gra A<r) = -p^(^ + ^), (5) 
Gpm p V dlnr d In r J 

where p is the mean molecular weight and m v is the proton 
mass. Numerical simulations (Evrard, Metzler & Navarro 
1996) show that this is a reliable mass estimator out to at 
least rsoo, the radius containing a mean overdensity of 500 
with respect to the critical density p c (z), which is compara- 
ble to the extent of emission r ex t in both groups (see below) . 
Taking p = 0.6 and using the emission-weighted mean tem- 
peratures given in Tabled along with the derived parame- 
ters for the gas profiles, the resulting total masses within r ex t 
are presented in Table|3] Adopting a Monte Carlo approach, 
the associated uncertainties were estimated by constructing 
10,000 artificial mass profiles. For each of these, the values of 
involved parameters (r c , /3, T) were taken from an asymmet- 
ric Gaussian distribution centered at the best-fitting value 
and having characteristic upper and lower widths given by 
the associated la uncertainties (Tables and |2J . The de- 
rived distribution of masses was in general also an asym- 
metric Gaussian; the quoted mass errors are the associated 
standard deviations (upper and lower) with respect to the 
nominal best-fitting values. 

Assuming the groups to be virialized out to T200, the 
virial radius can be self-consistently computed from the de- 
rived mass profile once p c {z) is known. Using 

po(z) = |^[(i + ^ 3 n m + (i + ^) 2 (i-n m -n A ) + fiA],(6) 

where the term in brackets is the ratio (H(z)/Ho) 2 , we de- 
rive nominal virial radii rw = r2oo = 0.57/i" 1 and 0.74/i _1 
Mpc for WJ943.7 and WJ943.5, respectively, in the adopted 
ACDM cosmology. We list r2oo and rsoo relative to r ex t in 
Table it is seen that emission has been detected beyond 
rsoo and out to 74 and 66 per cent of r2oo in the two cases. 

A temperature gradient in the X-ray gas changes total 
masses, and hence also r2oo and rsoo- Assuming the above 
polytropic gas profiles lowers the estimated masses within 
rsoo and r2oo by 9 and 20 per cent for WJ943.7 and 11 and 
29 per cent for WJ943.5 (rsoo and r2oo themselves decrease 
by 3 and 7 per cent for WJ943.7 and 4 and 11 per cent for 
WJ943.5). Given that the fractional errors on these masses 
would be larger than for the isothermal case (because of 
errors on 7), the two scenarios produce masses consistent 
with each other within la. Thus, within this statistical un- 
certainty the conclusions on M gra v remain unaffected by the 
presence of a temperature gradient. 
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4.3 Masses of X-ray gas 

To determine gas masses in the groups, the normalisation of 
the gas density profile must be established. For an isothermal 
ICM, this can be directly derived from the best-fitting 0- 
model for the surface brightness, using the spectral (in effect, 
luminosity) normalization in xspec, 

A = 10~ 14 (47r(D A (l + z)) 2 )- 1 J n e n H dV cm -5 , (7) 

where Da is the angular diameter distance to the source, 
and n e and nn are the number densities of electrons and 
hydrogen, respectively. With n e and uh following the best- 
fitting /3-profile out to r ext , their central values are found by 
evaluating the integral in eq. J7[, assuming n e /riH = 1.17. 
Wc find very similar central electron densities n e> o for the 
two groups of 3.15 ± 0.20 x 10~ 3 h 1/2 cm" 3 (WJ943.7) and 
2.91±0.19xl0~ 3 h 1/2 cm' 3 (WJ943.5), with errors obtained 
from the derived errors on the spectral normalization A. 
When converting these number densities to mass densities 
p — n e fi e m p — nH^HTrip, /i e = 1.17 and pn = 1-40 was 
assumed, appropriate for a fully ionized Z — 0.3^0 plasma 
(e.g. Mohr, Mathiesen & Evrard 1999). The total gas masses 
are finally found by simple volume integration of p(r) inside 

7°ext - 

In the non- isothermal case discussed above, M ga3 (r) 
changes because the derived gas density depends on emis- 
sivity and hence temperature. To quantify the magnitude of 
this change, we took for each group the best-fitting isother- 
mal spectral models within r ext and changed the associ- 
ated temperature within its extreme polytropic values at 
T(r = 0) and T(r = r ext ) (cf. Fig. g}, corrected to the 
'true' value via eq. @. This was done subject to the con- 
straint that the total 0.4-2.5 keV flux, from which the sur- 
face brightness profile has been extracted, should remain 
constant. The resulting change in spectral normalisation A 
then translates, via eq. (J7J, into a maximum resulting er- 
ror on the central gas density and hence total gas mass. In 
all cases, these changes were found to be less than 5 per 
cent within r ex t relative to the isothermal case, i.e. lying 
within the la statistical uncertainties from the spectral fits. 
We therefore assumed in the following that M gas remains 
unaffected by a change to the polytropic models. 



5 DISCUSSION 

5.1 Gas mass fractions 

The resulting gas mass fractions f gas = M gaa / M grav within 
r ext and rsoo are also presented in Table [3] along with 
the values resulting from extrapolating out to r2oo- The 
associated errors were obtained from M gas and M grav us- 
ing standard error propagation. For comparison with previ- 
ous work on the topic, we note that the derived gas mass 
fractions within (rsoo,?"20o) correspond to (0.049,0.065)/i _3,/2 
(WJ943.7) and (0.060,0.072)/i _3/2 (WJ943.5) for aflat stan- 
dard CDM universe, and (0. 053,0. 070)/i _3/2 (WJ943.7) and 
(0.067,0.079)ft" 3/2 (WJ943.5) for an open CDM universe 
with deceleration parameter qo — 0. Based on a coher- 
ently analysed sample of 66 systems ranging from individual 
galaxy haloes to massive clusters (T = 0.5 — 17 keV), Sander- 
son et al. (2003) derive a mean gas fraction within T2oo of 



<f g as>= 0.0785 ± 0.006h -3/2 (OCDM, q = 0), with, in 
general, lower values for cool systems and a statistically sig- 
nificant trend of f ga s increasing with T, a trend which is, 
however, subject to considerable scatter. The two groups 
studied here seem to conform to this trend, with (in this 
particular qo — cosmology) WJ943.5 lying close to this 
mean value and WJ943.7 having somewhat smaller f gas - In 
any given cosmology, the derived nominal gas mass fraction 
for the WJ943.5 group is also found to be comparable to 
mean values reported for more massive clusters (T > 4 keV) 
by Ettori & Fabian (1999) and Castillo-Morales & Schindler 
(2003) and for intermediate-mass systems (T = 2.5 — 5 
keV) by Mohr et al. (1999) and Arnaud & Evrard (1999). 
WJ943.7, however, displays a lower f gas than found in these 
studies, with a value very similar to, e.g., that found for the 
T ~ 2 keV cluster A1983 (Pratt & Arnaud 2003). For both 
groups, f gas is found to be clearly smaller than the values 
obtained for massive (T > 4 keV) clusters (e.g., Mohr et 
al. 1999; Arnaud & Evrard 1999; Grego et al. 2001; Sander- 
son et al. 2003). In summary, within r2oo WJ943.5 seems 
to be reasonably representative of the general cluster pop- 
ulation, whereas f ga3 of WJ943.7 is more typical of values 
reported for slightly less massive systems situated at the 
interface between groups and clusters. 

Using the polytropic models, we find that f ga3 is raised 
to 0.061 (0.086) within r 50 o (r 2 oo) for WJ943.7, and 0.079 
(0.107) for WJ943.5. These values are larger than the nomi- 
nal isothermal values by 5 per cent within rsoo and between 
10 and 20 per cent within r2oo- Even more so than for the 
total mass M grav , they are comfortably consistent with the 
isothermal values within la. We stress that, barring any 
additional systematic errors, our derived values within r2oo 
should be fairly reliable, since they involve a relatively small 
degree of radial extrapolation compared to results for other 
groups. In this light it is interesting to note that our val- 
ues are much closer to cluster values than often reported for 
small groups within smaller radii (e.g. David, Jones & For- 
man 1995; Pedersen, Yoshii & Sommer-Larsen 1997; Mohr 
et al. 1999; Mushotzky et al. 2003). 

The modelled radial variation of f gas in both groups is 
shown in Fig. |H] For the polytropic models, M grav is lower 
at large radii, but rsoo and r2oo, within which M gas is eval- 
uated, also decrease. Since f gaa clearly increases with r even 
in the isothermal case, the relative impact of assuming a 
negative temperature gradient is therefore less pronounced 
at large radii for f gas than for M grav The gas mass frac- 
tion of both systems increases out to the radii probed, with 
the standard implication that the intragroup gas is more ex- 
tensively distributed than the dark matter. The consistency 
between the isothermal and polytropic models strongly sug- 
gests that when taking the errors into account, the overall ra- 
dial behaviour of f gas shown in Fig.|S|is a robust result. The 
strong dependence of f gas on radius suggested by the figure, 
particularly for the cooler WJ943.7, is furthermore corrob- 
orated by the results of *l3.2l for isothermal NFW fits to the 
surface brightness profiles, which, when scaled to r c or r 2 oo, 
show NFW scale radii for the gas much larger than those 
of simulated and observed hotter systems with a milder ra- 
dial variation of f gaa (cf. Ettori & Fabian 1999). The weak 
temperature dependence of the emissivity in the relevant 
temperature and energy ranges ( q4.1H implies that, within 
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Table 3. Derived gas masses, total gravitating masses, and gas mass fractions within the derived extent of emission r ex t, rsoo, and the 
approximate virial radius r200- All values are given in the adopted ACDM cosmology of (Q m ,fl/^) = (0.3,0.7). 
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Figure 6. Gas mass fraction f ga s plotted from r c to r200, us- 
ing two assumptions for the temperature profile. Dark curves for 
WJ943.7, light curves for WJ943.5. Solid lines represent isother- 
mal gas, dashed lines the polytropic models discussed in the text. 
Typical error bars for the isothermal case are shown at r = r e xt, 
beyond which fgas is extrapolated. 



the errors, this result holds equally well for the polytropic 
models. 



5.2 Entropy 

The entropy of the ICM, here defined as S = T jn'J^ , reflects 
the accretion history of the hot gas and as such is a prime 
indicator of the physical processes involved in the formation 
and evolution of groups and clusters. Recent studies of the 
entropy in systems covering a wide range in masses have 
made it clear that gravity-driven processes alone cannot ac- 
count for the observed entropy: Systems at all scales, most 
notably low-mass ones, show evidence for excess entropy rel- 
ative to that expected from purely gravitational heating and 
compression of the gas. 

Various mechanisms have been proposed to explain this 
excess entropy, see e.g. the discussion in Ponman, Sanderson 
& Finoguenov (2003). We note here that the rather large 
values of f gas found for our systems seem inconsistent with 
cooling-only models, which require more than 50 per cent of 
the gas in low-T systems like these to cool out of the diffuse 
hot phase (e.g. Dave, Katz & Weinberg 2002; Ponman et al. 
2003). 

The radial variation of entropy cannot be computed in 



a model-independent way for the two groups studied here, 
given the poorly resolved temperature profiles. Using instead 
the best-fitting /3-models along with the isothermal assump- 
tion, resulting central entropies are 78 ±8 and 118 ±20 /i~ 1//3 
keV cm 2 for WJ943.7 and WJ943.5, respectively, rising to 
113 and 142 h~ 1/3 keV cm 2 at a fiducial radius of 0.1r2oo. 

Self-similar scaling between systems of different masses 
would lead to S ex (1 + jz) _2 T. In Fig.[7| modelled profiles of 
S(r) are plotted for each group for both the isothermal and 
polytropic cases. The entropy has been scaled by (1 + z) 2 /T 
so as to 'remove' the dependence on mass and redshift in 
the self-similar case. In this representation the profiles of 
the two groups are seen to be fairly similar. The flattening 
in the entropy profiles seen inside 0.1r2oo, due to the core in 
the gas density, is not a reliable result, since it is based on 
the assumption that T(r) does not drop within the central 
~ 50 kpc. As can be seen from Fig. [I] our data do not allow 
us to resolve the temperature structure on such small scales. 
What is clear from our results, is that any isentropic core 
in these groups does not extend to r>0.2r2oo, in contrast 
to the predictions of simple preheating models for such poor 
systems (Tozzi & Norman 2001, Babul et al. 2002; Voit et al. 
2003). This is consistent with recent results for other groups 
(Mushotzky et al. 2003; Pratt & Arnaud 2003) and argues 
against significant preheating in a simple spherical collapse 
scenario. 

Also plotted in Fig. |7| (triangles) is the expectation for 
S based on self-similar scaling of massive clusters (Ponman 
et al. 2003). For both groups, the scaled entropy exceeds 
that observed in massive clusters at 0.1r2oo- Further out 
our results are more uncertain, but the modelled profiles 
suggest that this excess seems to persist out to at least rsoo, 
where the scaled (isothermal) entropies of the groups are a 
factor ~ 2 above the self-similar values. For our observed gas 
density profiles, the polytropic indices of both groups would 
have to be 7 > 1.40 in order to avoid excess entropy at rsoo- 
This is well beyond the values of 7 obtained for observed and 
simulated groups and clusters. The derived entropy levels 
at 0.1r2oo and rsoo, and hence the associated excesses with 
respect to the massive clusters of Ponman et al. (2003), are 
in good agreement with the best-fitting S — T relations at the 
relevant radii of the entire group/cluster sample of Ponman 
et al. (2003). 

Spherical shock heating models and most cosmologi- 
cal simulations produce entropy profiles with a power law 
form S oc r (e.g. Tozzi & Norman 2001). Two lines of 
this slope have also been plotted in Fig. |7| The lowest line, 
which agrees fairly well with the triangles representing ob- 
servations of rich clusters, is normalised to the scaled mean 
value at r2oo of the M ~ 5 x 10 14 hT^M® simulated 'Santa 
Barbara' cluster (Frenk et al. 1999), which included only 
gravitational physics. The profile has been scaled using an 
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Figure 7. Scaled entropy S(l + z) 2 /T as a function of r/r200, 
using two assumptions for the temperature profile. Triangles mark 
the expected values at 0.1r2oo and r 500 based on a self-similar 
scaling from the values observed in hot clusters. Sloping dotted 
lines represent the self-similar profile expected from shock heating 
in a spherical smooth accretion scenario, S oc r , normalised at 
f*200 to the preheated T fs 2 keV cluster simulation of Borgani 
et al. (2001), and to the simulated T pa 6.5 keV cluster of Frenk 
et al. (1999) which includes no heating or cooling effects. 



estimated emission-weighted temperature of T ew — 6.5 keV, 
derived from its mass-weighted value of T mw — 4.7 keV as- 
suming a ratio T ew /T em ~ 1.4, as suggested by results from 
simulations (Borgani et al. 2002) and observations (Sander- 
son et al. 2003). While only a suggestive result, this pro- 
file clearly fails to reproduce our modelled entropy levels at 
large radii. Also shown is an S oc r 1 ' line normalised at 
T200 to the result obtained for a T = 2.1 keV cluster in nu- 
merical simulations of cluster formation by Borgani et al. 
(2001). These simulations include preheating at a level of 
45/T 1/3 keV cm , and do reproduce the large entropy ex- 
cesses hinted at for our groups at rsoo without forming a 
distinct isentropic core. 

Whereas the absence of a large isentropic core is in di- 
rect conflict with simple smooth accretion of preheated gas, 
the excess entropy at large radii suggested by our analyti- 
cal profiles may, in fact, be explained by such models (for 
details, see Voit et al. 2003 and Voit & Ponman 2003). 



6 CONCLUSIONS 

The main results of this work may be summarized as follows. 

• Using sensitive XMM-Newton observations, we have 
mapped the X-ray emission out to, on average, 70 per cent 
of the estimated virial radii V2oo in two X-ray selected galaxy 
groups at redshifts z = 0.18 and z — 0.256. We find the X- 
ray surface brightness profiles of both systems to be well de- 
scribed by standard /3-models out to the radii probed, with 
indices of (3 — 0.49 and 0.63. Fits to isothermal NFW profiles 
for the gas density indicate that the surface brightness data 



do not rule out a steepening in the density profile at large 
radii as might be expected on physical grounds. However, 
adequate fits for such models can only be obtained with a 
very large scale radius, and models in which gas traces mass 
are strongly ruled out for both groups. 

• Emission-weighted mean temperatures, determined 
within the same radii as the surface brightness profiles, are 
found to be kT = 1.7 ± 0.1 and 2.4 ± 0.4 keV, with cor- 
responding gas metallicities roughly 0.3 solar. The temper- 
ature profiles are poorly resolved but are consistent with 
isothermality within 90 per cent confidence. 

• Radial profiles of gas mass, gravitating mass, and gas 
mass fraction have been determined assuming both isother- 
mal and polytropic gas distributions allowed by the data. In 
both cases, the gas mass fraction clearly increases with ra- 
dius, implying a more extended gas distribution than that of 
the total mass. Within the derived extent of emission, gas 
mass fractions are 0.065h~ 3 ^ 2 and 0.077/i -3 ^ 2 in the two 
groups, considerably higher than values reported in the lit- 
erature for other groups within smaller radii. Extrapolated 
gas mass fractions at r2oo of 0.077fo _3//2 and 0.089fa -3//2 are 
slightly lower, but comparable to, values found for more 
massive clusters. These results emphasize the importance 
of extending group studies to large radii, in order to ob- 
tain reliable information on the amount and distribution of 
baryonic and dark matter in these systems. 

• The radial variation of gas entropy in the two groups 
shows no evidence for a large isentropic core, at vari- 
ance with predictions from simple preheating scenarios. The 
poorly resolved temperature profiles prohibit any firm con- 
clusions on the entropy behaviour at large radii, but evi- 
dence for excess entropy with respect to that expected from 
self-similar scaling of observed massive clusters is seen out 
to at least rsoo- Such a large-radius excess would be consis- 
tent with the predictions of some simulations including heat 
sources, and with smooth accretion theory. As an alterna- 
tive for the origin of excess entropy at all scales, radiative 
cooling alone is ruled out by the relatively high gas mass 
fractions derived for these systems. 
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